Probing nucleon-nucleon interactions in breakup of a one-neutron halo n Be on a 

proton target 
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A comparison between full few-body Faddeev/Alt-Grassberger-Sandhas (Faddeev/AGS) and 
continuum-discretized coupled channels (CDCC) calculations is made for the resonant and non- 
resonant breakup of n Be on proton target at 63.7 MeV/u incident energy. A simplified two-body 
model is used for n Be which involves an inert 10 Be(0 + ) core and a valence neutron. The sensitivity 
of the calculated observables to the nucleon-nucleon potential dynamical input is analysed. We 
show that with the present NN and N-core dynamics the results remain a puzzle for the few-body 
problem of scattering from light exotic halo nuclei. 
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Halo nuclei are few-body light nuclear systems that ap- 
pear at the neutron drip line. The study of these nuclei 
is providing new theoretical challenges to the nuclear re- 
action theory, a key tool to interpret and extract reliable 
nuclear structure information from the new generation of 
very precise experimental measurements. 

When describing the scattering of stable nuclei from 
halo nuclei it is crucial to handle its few-body character. 
In addition it is necessary to treat all opening channels 
(elastic, inelastic, transfer and breakup) on equal footing. 
Yet, it is fair to say that a tight control on the under- 
lying physics of the reaction mechanisms has not been 
achieved. 

We aim to shed light on the relevant dynamical aspects 
of the reaction framework. 

Recently the resonant and non-resonant breakup re- 
sulting from the scattering of 11 Be from a proton tar- 
get at 63.7 MeV/u was measured at MSU in order to 
obtain information on the 11 Be continuum by nuclear 
excitation [![. The experimental results were compared 
with the calculated observables using the Faddeev/Alt- 
Grassberger-Sandhas (Faddeev/AGS) scattering frame- 
work [2j-[4j and making use of a simplified two-body model 
for Be in terms of an inert 10 Be(0 + ) core and a valence 
neutron Due to the experimental energy resolution, 
the breakup angular distribution results contained inte- 
grated contributions from ranges of relative neutron-core 
energy. It was found in Q that in the case where the rel- 
ative core-neutron energy is integrated around the reso- 
nance E T = 1.275 MeV in the energy range E IC \ = — 2.5 
MeV the Faddeev/AGS approach reproduced fairly well 
the shape distribution of the data although underesti- 
mating the maximum value of the breakup observable. 
A large discrepancy was found between the predictions 
of the Faddeev/AGS approach and those made by the 
continuum-discretized coupled channels (CDCC) frame- 
work done in the works of Q], [||. Given the previous 



p+ 11 Be benchmark calculations [7J one should expect 
some differences between Faddeev/AGS and CDCC re- 
sults, however, the different dynamic input may be re- 
sponsible for the disagreement between the predictions 
of (| and [J, i] as well. In this work we aim to make a 
fair comparison between the two scattering frameworks 
using the same dynamical and structure inputs in or- 
der to establish the source of the discrepancy and shed 
light on the sensitivity of the observables to some key 
aspects of the dynamic input. In any 3-body reaction 
approach to our working example of the breakup of 11 Be 
by a proton target, one needs as input the 3-pair poten- 
tials: nucleon-nucleon (NN), N-core and p-core. In the 
late nineties high precision NN potentials wich reproduce 
the NN data below 350 MeV lab energy with \ 2 /datum 
~ 1 were developed such as the CD Bonn Q and AV18 
potentials. However, puzzling 3 body problems re- 
main (such as the Ay discrepancy) where the current NN 
potentials are unable to reproduce the data. It is relevant 
to know to what extent our three-body observables are 
sensitive to this underlying dynamical input. 

The Faddeev/AGS is a non-relativistic three-body 
multiple scattering framework that can be used to cal- 
culate all relevant three-body observables. We use here 
the odd man out notation for the three interacting parti- 
cles (1,2,3) which means, for example, that the potential 
between the pair (2,3) is denoted as v\. According to this 
reaction framework, one needs to evaluate the operators 
U 13 ", whose on-shell matrix elements are the transition 
amplitudes. These operators are obtained by solving the 
three-body AGS integral equation [H, 0| 
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with a, (3, 7 = (1, 2, 3), (j3 — in the final breakup state). 
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Here, S/s a = 1 — 5p a and the pair transition operator is 

i 7 = v 1 + VjGoty , (2) 

where Go is the free resolvent Go = (E + iQ — Hq) -1 , and 
E is the total energy of the three-particle system in the 
center of mass (cm.) frame. 

In our calculations Eq. (TTJ) is solved exactly in mo- 
mentum space after partial wave decomposition and 
discretization of all momentum variables. The Pade 
method [lOj is used to sum the multiple scattering se- 
ries. We include the nuclear interaction between all three 
pairs, and the Coulomb interaction between the proton 
and 10 Be, follo wing the technical developments imple- 
mented in Refs. |lll.ll2j| and the breakup observables are 
calculated as summarized in detail in [13| . 

The 3-body Continuum-Discretized Coupled-Channels 
(CDCC) \l4 . [HI reaction framework consists in solv- 
ing the Schrodinger equation in a model space in which 
the three-body wavefunction is expanded in the internal 
states (bound and continuum resonant and non-resonant 
states) of the two-body projectile H p . In practical cal- 
culations, the continuum spectrum has to be truncated 
in excitation energy, and discretized into a set of square- 
integrable functions. The most widely used discretization 
method is the called binning method, in which the con- 
tinuum is divided into a set of energy intervals; for each 
interval, or bin, a representative function, (p a (r), is con- 
structed by superposition of the true scattering states 
within the bin interval. The total three-body wavefunc- 
tion is expanded in terms of these representative func- 
tions as 

N 

^ cc (R,r) = 5> a (r)u; a (R) (3) 

where in here a=0 refers to the projectile ground state, 
Ko is the incident wave number of the projectile in the 
cm. frame, R the relative distance between the cm. of 
the projectile and the target and r the relative distance 
between the valence particle and the core. The wave 
functions uj a (R) of the projectile-target relative motion 
are solutions of the coupled-channel equations 

(E a -T R - V aa (R)) w a (R) = ^(R)^(R) , (4) 

where E a — E — e a and the coupling potentials are 

V aP (R) = (<j> a \ ]T ^(R,r)|^> . (5) 

j=C,v 

The CDCC method was originally developed in order to 
incorporate the effect of the breakup channels in deuteron 
induced reactions. The proton-target and neutron-target 
interactions used in these calculations are usually taken 
as optical potentials adjusted to reproduce the elastic 
scattering at the same energy per nucleon, i.e., E p « 
E n k, Ed/2. For nucleoli- nucleus scattering these optical 



potentials are in many cases well represented by a sim- 
ple, local, L-independent interaction comprising a cen- 
tral and, maybe, a spin-orbit term [l7|. This has been 
also the standard choice in the application of the CDCC 
method to the scattering of other weakly bound nuclei 
( 6)7 Li, n Be, 8 B, etc) by medium -mass or heavy targets. 
However, this simple prescription might not be appropri- 
ate to describe the scattering of halo nuclei on protons, 
because in this case one of the fragment-target interac- 
tions is the NN potential, which is known to be strongly 
L-dependent. Furthermore, the absence of an imaginary 
(absorptive) part in the NN interaction makes less clear 
the formal justification of the CDCC method as an ac- 
curate approximation of a three-body scattering problem 

In this context, the calculations presented in this work 
arise from a twofold motivation. Firstly, we aim to study 
the importance of using a realistic NN interaction in the 
description of the scattering of one-neutron halo nuclei by 
a proton target. In addition, by comparing the CDCC 
with the Faddeev/AGS calculations using the same three- 
body Hamiltonian we check to what extent the CDCC 
method provides an accurate reaction framework to the 
solution of a three-body scattering problem, following 
previous work done in ref. Q. 

In order to study the sensitivity of the calculated ob- 
servables to the underlying NN potential, we have com- 
pared the Faddeev/AGS calculations using the realistic 
CD Bonn d| and AV18 [9] NN potentials with those ob- 
tained making use of a simple L-independent potential. 
For the latter, we first consider a simple Gaussian poten- 
tial v pn (r) = -w exp[-(r/r ) 2 ], with v Q = 72.15 MeV 
and ro = 1.484 fm. These parameters are adjusted 
to fit the deuteron binding energy and low-energy 3 Si 
proton-neutron phase-shifts. This parametrization, here 
denoted as Gte, has been used in several works to gener- 
ate the deuteron states in d+A CDCC calculations (see 
e.g. Ref. [IBl). The scattering phase-shifts associated to 
this potential for S- and P- waves are compared in Fig. Q] 
with those obtained with the realistic CD Bonn potential. 
It is seen that the 3 Si phase-shifts are well described by 
the Gaussian parametrization up to a cm. energy of » 20 
MeV, but they are strikingly different from the realistic 
phase-shifts for the singlet S-wave and the P-waves. 

In addition, we have performed some additional Fad- 
deev/AGS calculations considering several L-dependent 
modified Gaussian potentials: the Gse (Gpi) where the 
S (P) partial waves were modified taking the parameters 
from a Gaussian potential adjusted to reproduce the sin- 
glet scattering length [l6[ and all the other partial waves 
kept the same as the Gte potential. Finally, we also con- 
sider the Gp2 example case where the P and P-F partial 
waves were taken from the realistic CD Bonn potential 
and all the other partial waves kept the same as the Gte 
potential. The proton-core and neutron-core pair inter- 
actions are taken as described in Ref. [|[ . 

In the solution of the Faddeev/AGS equations we in- 
clude n-p partial waves with relative orbital angular mo- 
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FIG. 1: (Color online) NN phase shifts for S and P partial waves obtained from the CDBonn potential (dashed lines) and 
Gaussian potential (solid lines) as indicated by the labels, where the depth and range parameters were adjusted to fit the 
deuteron binding energy, and denoted in the text as Gte (light lines) 
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FIG. 2: (Color online) Energy spectrum for p( 11 Be,p) 10 Be n at 63.7 MeV/u integrated over the whole angular range. The NN 
interactions are described in the text. 



mentum L < 6, p- 10 Be with L < 21, and n- 10 Be with 
L < 6. Three-body total angular momentum is included 
up to 40. 

As for the CDCC equations we discretized the n- 10 Be 
continuum using energy bins for partial waves L < 4, and 
up to a maximum excitation energy of 40 MeV, above 
breakup threshold. The coupled equations were solved 
for total angular momentum up to 60, and the solutions 
matched to their asymptotic form at a distance of 60 
fm. Both Coulomb and nuclear couplings were included 
and treated on equal footing. For the energy spectrum 
the number of bins was increased for excitation energies 
below 2 MeV in order to provide a finer description of 
the resonance region. 

In Fig. [5] we represent the calculated energy spectrum 
da/dE le i that emerges by integrating the semi- inclusive 
cross section over the solid angle dO c m- . As follows 
from Fig. f2^) the Gaussian L- independent potential Gte 
(dashed line) considerably overpredicts the non-resonant 
background and also the resonant peak with respect to 
the Faddeev/ AGS result calculated with the CD Bonn 



potential (solid line). The predicted results using CD 
Bonn are indistinguishable from those which make use 
of the AV18 @ potential and are not represented in the 
graph. The CDCC result (dashed-dotted line) that uses 
the same Gt e potential is very close to the corresponding 
Faddeev/ AGS calculation except at the very low relative 
energies in a region around 1 MeV, where the CDCC 
cross section is somewhat smaller. The sensitivity to the 
underlying NN low partial waves interactions is shown in 
Fig. fSb) which shows that this observable is sensitive to 
a relistic treatment of the NN potential in particular of 
the P- waves. 

In Fig. [3] we show the breakup angular distribution 
d<j/d£l c . m .. We have not included very small angles 
^c.m. < 5° where there are no data and the convergence 
of the Faddeev/ AGS results with respect to the Coulomb 
screening radius is slow. Due to the energy resolution of 
the experimental setup, the relative core- neutron energy 
is integrated around the resonance E r — 1.275 MeV in the 
energy range E ie i = — 2.5 MeV in the upper part of the 
figure. The angular distribution calculated with the Fad- 
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FIG. 3: (Color online) Angular distribution for the breakup p( 11 Be,p) 10 Be n at 63.7 MeV/u integrated over the energy range 
E Te \ = — 2.5 MeV (upper part) and integrated over the energy range E Te \ — 2.5 — 5.0 MeV (lower part). The NN interactions 
are described in the text. 



deev/AGS approach and using the CD Bonn NN poten- 
tial (solid curve) shown in Fig. [3^) describes the overall 
shape distribution of the data reasonably well, although 
clearly underestimating the data around C m = 20 de- 
grees by about 40%. The dependence on the calculated 
observable on the proton-core optical potential was stud- 
ied. Using a potential that fits the elastic proton core 
data 0| increases the calculated observable at 8 c . m . = 20 
degrees by less than 10%. The origin of this disagreement 
needs to be further investigated. The Faddeev/AGS re- 
sults are again indistinguishable when using the AV18 
Q potential. It follows that these 3 body observables are 
insensitive to the choice of the realistic NN potential and 
probe essentially the NN scattering on-shell behaviour. 

Upon substituting the realistic potential by a L- 
independent Gaussian Gte potential the Faddeev/AGS 
results for the calculated angular distribution (dashed 
solid line) are significantly enhanced at small cm. an- 
gles and overall become similar to the predictions of 
the CDCC approach (dashed-dotted line). As seen in 
Fig. [SJs), in this case, the angular distribution is very 
sensitive to a realistic treatment of the NN potential, in 
particular to the NN P-waves. 

In the lower part of Fig. [3] we show the breakup angular 
distribution da / 'rfO c m . , where the relative core- neutron 
energy is integrated over the energy range E Te \ = 2.5 — 5.0 
MeV. In this case the Faddeev approach using the CD 
Bonn potential does not reproduce the data. As in 



the upper case the Faddeev calculation using the GTE 
Gaussian interaction overestimates the realistic calcula- 
tion and overall becomes similar do the CDCC approach. 
Again the nonresonant breakup at higher relative ener- 
gies is also very sensitive to a realistic treatment of the 
NN low partial waves. 

Overall we can say that in our case study, there is a 
fairly good agreement between the two microscopic reac- 
tion formalisms, the Faddeev/AGS and the CDCC, bet- 
ter than in the benchmark calculation of ref. (tJ. We 
attribute this difference to the fact that, in the referred 
work, the breakup observables are studied with respect to 
the 10 Be core (integrated with respect to the neutron an- 
gle). That calculated inclusive breakup observable, was 
found to be dominated by configurations with small p-n 
angular momentum. These configurations are difficult to 
treat in a CDCC framework based in the expansion of the 
internal states of the n- 10 Be [Eq. [3J, and hence the slow 
convergence found in Ref. [7] for those observables. In 
the present work, by contrast, we study selected regions 
dominated by small n-core relative energies and angular 
momenta, for which the convergence of the CDCC ap- 
proach is expected to be much faster. Similar conclusions 
were achieved in the work of Ref. [19j for the breakup of 
8 B on the 58 Ni target. 

In conclusion, we have performed full Faddeev-type 
and CDCC calculations for the breakup of 11 Be on a pro- 
ton target at 63.7 MeV/u incident energy. A simplified 
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two-body model for n Be consisting of an inert lo Be(0 + ) 
core and a valence neutron has been used. We have 
shown that the Faddeev results using a i-independent p- 
n potential of Gaussian form (with the depth and range 
parameters adjusted to fit the deuteron binding energy) 
considerably differ from those obtained with a realistic 
CD Bonn potential. The former are found to overes- 
timate the resonant and non-resonant breakup angular 
distributions calculated with the realistic NN interaction. 
These results strongly suggest that, at least in this en- 
ergy regime, these breakup observables are very sensitive 
to the underlying NN interaction, particularly to the P- 
wave. 

We have also compared the Faddeev calculations with 
CDCC calculations, using in both cases the simple NN L- 
independent Gaussian parametrization. The angular and 
energy breakup distributions are fairly similar, display- 
ing only some small differences in the magnitude of the 
breakup cross section at small excitation energies. The 
good agreement between the two formalisms in this case 
study is better than in the previous benchmark calcula- 
tion Ref. 0. 

Since existing CDCC codes usually rely on simple cen- 



tral, ^-independent fragment-target interactions, exten- 
sions of these codes, in order to incorporate realistic NN 
potentials in the calculation of the coupling potentials 
would be mandatory for future applications of the CDCC 
method in order to extract physically meaningful infor- 
mation from nucleon-nucleus scattering data. 

We conclude that one can only extract reliable infor- 
mation from the breakup of halo interaction if a realistic 
potential is used. Different realistic NN potentials pre- 
dict the same breakup observables which then essentially 
probe the NN scattering on-shell behaviour. In addition, 
the breakup angular distribution integrated around the 
resonance underestimates significantly the data by about 
40% and clearly does not reproduce this observable when 
integration is made at a higher relative core-neutron en- 
ergy range. Further work should be performed to under- 
stand the physics of this discrepancy. With the present 
NN and N-core dynamics these results remain a puzzle 
for the few-body problem of scattering from light exotic 
halo nuclei. 
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